Statistical thermodynamics and weighted topology of radial networks 
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We propose a method of characterizing radial networks based on a partition function associated 
with the structural triangulation of the network. The internal energy, Helmholtz free energy, and 
entropy derived from the partition function are used to group similar networks together and to 
interrogate the history of their development. We illustrate our methodology for a model for optimal 
transport in tree leaves. 
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Interest in networks has expanded rapidly in recent 
years [T]-[5]. In particular, there has been substantial 
interest in problems associated virith modelling venation 
patterns [IHS] , optimizing supply networks [71 |H] , and de- 
termining final system shapes and sizes based on a sim- 
ulated growth mechanisms ^ . 

The characteristics of a network are determined by the 
supply and demand mechanisms of the system. For ex- 
ample, a leaf is supplied with nutrients and sunlight as 
energy sources, and demand is created by the biological 
drive to keep the leaf cells alive and the leaf growing. For 
a river ending in a delta, the supply is the amount of rain 
falling, and demand is created by the pull of gravity, the 
amount of water dissipation, and the terrain of the delta. 

The rules and parameters governing the supply and 
demand of a system are expected to be invariant over 
time scales beyond its lifespan. Hence, the same rules and 
parameters govern the creation and the extinction states 
of a system. What happens between these is affected by 
external factors. For example, a leaf may suffer under less 
than ideal weather conditions or a river may be affected 
by the construction of a dam. 

An important question is whether systems can be 
grouped according to their shape and venation at a given 
state in their lifespan. Going a step further, one may ask 
whether interpretations of life history can be made by 
analyzing and comparing a state in the lifespan of a sys- 
tem to a comparable state of another system. We expect 
that the growth history of a system is encoded in the dis- 
tribution of energy in its venation, i.e. the work done in 
building up the supply system. Due to external factors 
and the existing structure of the system, some parts of a 
venation system may have required more energy to build 
up than other parts. 

In this paper, we introduce a method for measuring the 
energy distribution of the venation of a network based 
on a partition function calculated from its energy states, 
as represented by a structural triangulation. The par- 
tition function is used to calculate thermodynamic vari- 
ables, which enable systems to be grouped together at 
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their maturation state - the state when the system stops 
growing, as well as enabling the development of a net- 
work to be investigated. We illustrate our methodology 
by analyzing simulated tree leaves. 

Topological invariants are used to characterize shapes 
and structures of geometrical objects. A prominent ex- 
ample is the Euler characteristic x |10j - a measure of 
the topological structure of a space that is invariant un- 
der deformations. For convex polytopes, for example, 
X^V — E + F = 2 where V, E and F are the numbers 
of vertices, edges and faces, respectively, of the polytope. 

We will focus on planar radial networks with a sin- 
gle origin and without loops. Let the venation of the 
system consist of edges e = {vi,Vj) connecting vertices 
Vi = {xi,yi), with one vertex designated as the origin 
Wo = (0,0), which is associated with the first generation 
of the network, g{vo) = 1. Then, an edge e = {vi,Vj) 
always connects an nth generation vertex, g{vi) = n, to 
an (n -I- l)th generation vertex, g{vj) = n + 1. 

We define a triangular element A as a set of three 
vertices, A = {vi^Vj,Vk}, such that the generational as- 
sociation is g{vi) = n and g(vj) = g{vk) = n + 1. A tri- 
angular element is assigned the generation <?(A), which 
corresponds to the lowest generational number of its ver- 
tices. We emphasize that there is no intersection between 
triangular elements. The set of all triangular elements in 
a network is called the triangulation A4 of that network. 
Although we have not proven so, the uniqueness of the 
triangulation of a network according to the foregoing pre- 
scription is a plausible working hypothesis. 

Two venation networks and their triangulations are 
shown in Fig. [T] Each triangular element A^ not only 
has an associated generation g{Ai), but also an associ- 
ated area A{Ai). While the structure of a triangulation 
is unchanged by deformations of the network, the areas 
of the elements allow the original and transformed net- 
works to be distinguished. Accordingly, we use the gen- 
eration and area associated with the triangular elements 
as a measure of the energy distribution in a venation. In 
effect, we have a weighted topological description of the 
network based on its triangulation. 

For increasing generation number g, the system ex- 
pends correspondingly more energy creating a vertex v 
and element A associated with g. The system must move 
material from the supply source at the origin to the po- 
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FIG. 1: (Color online) Two venation networks (upper panels) 
and their triangulations (lower panels). 



sition of g. A reasonable physical interpretation for g 
is, therefore, an energy level of the system. If the en- 
ergy unit of the system is -Eqj the energy Eg is the work 
required to produce a unit at generation g: Eg — gEg. 
Building on this point of view, a natural interpretation 
of the area A{A) of an element A is a measure of the 
material needed to assemble this element at generation 
f/(A). In general, if the material needed for an element 
of the system is N, then N{A) = pA{A), where p is a 
proportionality constant. Without loss of generality, we 
take A^ to be a positive integer, and set p = 1- This re- 
quires the coordinates of the vertices to be integer- valued 
as well, Xi S Z, such that the network is embedded in an 
integer lattice Z^. This modification is not a requirement, 
but it allows us to interpret N{A) as the number of 'par- 
ticles' with energy Egi^i^) required to create the triangle 
element A. 

We are now able to define the partition function of the 
triangulation A^ as 



Z{M) = Y^ N{A) ts 
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where t = e~^^° and /3 the inverse temperature. The 
factor t^^^^ originates from the fact that particles in the 
same generation are assigned to the same energy level 
-^g(A) = 5(^)-£'0- The partition functions for the net- 
works in Fig. [T] are 

Z{Mc) = 16t 



9*2 + Qt^ + 



4*4 + ^5 , 



(2) 
(3) 



where Aic and Md are the triangulations in Fig. IIJc) 
andllrd), respectively. 

The partition function is the sum over the microstates 
the system can occupy. These are identified in the net- 
work as 'particles' that have an associated history. To 
create a particle in element A, the network had to oc- 
cupy a specific state where work £'g(A) was done to create 
this particle. Accordingly, when we sum over all network 
system particles in (fTl), we are indeed summing over all 
microstates in its history with associated energy -Eg(A)- 
Because energy levels depend on the generation g of the 
triangle elements which, in turn, depend on the structure 
of the venation, the partition function in (fl]) can be con- 
sidered not only as a statistical measure of the network, 
but also as a measure of the structure of the venation. 
Therefore, we expect two systems with similar partition 
functions to display similar venation structures. 

The partition function enables us to calculate thermo- 
dynamic variables for the network. The internal energy 
is 



U 



d\nZ 
dp 



(4) 



Figure[2]^a) shows that U{Mb) > U{Ma) for all /3. In the 
first few generations (large /3), the internal energies of the 
two networks are similar, but diverge as higher genera- 
tions are included (decreasing j3). The physical interpre- 
tation of this is that there has been more work done by 
the network system to create the venation of the asym- 
metric network in Fig. m[b) than that of the symmetric 
network in Fig. fl|a). Thus, the structure of the asym- 
metric network displays a higher internal energy through 
more elements being in higher generations of the system 
than in the symmetric network. These qualitative trends 
can already be seen in the partition functions. In general, 
we expect that a network which displays a higher degree 
of venation symmetry than another network system to 
have the lower internal energy. 

The entropy S — fc^ (In Z -I- f5U) in the canonical en- 
semble, where we have set Boltzmann's constant fcs = 1. 
As Fig. [gb) shows, S{Ma) > S{Mt) for all /3. To un- 
derstand this result, we identify entropy with the energy 
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FIG. 2: (a) Internal energy and (b) entropy for the networks 
in Fig. ITJa) (solid curves) and Fig. [T|b) (broken curves). 




FIG. 3: (Color online) Simulated leaves (top panels) with parameters (a) a — 0.6, j3 = 0.5, (b) a — 0.65, jS = 0.38, (c) 
a = 0.6, /3 = 0.8, (d) a = 0.75, P = 0.7, and (e) a = 0.65, /3 = 0.7, with corresponding triangulations (bottom panels). All 
leaves have e — 4.0. 



dispersal in the network, rather than as a measure of 
structural 'disorder' of the venation. The energy in the 
symmetric network is more evenly distributed between 
the different generations than in the asymmetric network 
and, hence, has the higher entropy. Two limiting forms 
of the entropy are especially useful. As /3 — >■ oo (i — > 0), 
the entropy approaches the logarithm of the triangulated 
area associated with the first generation: 



lim 5' = In 



N{A) 



9=1 



(5) 



The other limit is /3 -^ (i ^ 1), where the entropy 
approaches the logarithm of the area of the entire trian- 
gulation: 



lim S 
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Both of these limits can be obtained directly from the 
partition function; the interesting situation is when the 
entropy curves of different networks cross. 

We consider now a model for leaf growth based on 
the notion of optimized transport. Supplying a leaf with 
water and nutrients has a cost. Minimizing that cost 
determines the shape of the leaf. We follow Xia 9J in 
developing a model for the supply network of tree leaves 
that is parametrized by three quantities that character- 
ize the cost function associated with extending the net- 
work. The transport cost is minimized when nutrients 
are moved in the same direction and maximized when 
moved in opposing directions. We eliminate the option 
of opposite flows in the transport network by setting the 
cost of sustaining such a flow to infinity. For all other 



flows, the parameter /3 scales the extent to which the flow 
direction contributes to the total cost of sustaining the 
supply network. Another cost is associated with supply 
flowing through a network edge. The more matter that 
flows through an edge, the larger the cross-sectional area 
of that edge which requires more energy to sustain the 
pressure in the descendant supply network. The 'thick- 
ness' of a given edge is proportional to its weight, so 
the cost of maintaining an edge with a particular weight 
leads to a second parameter, a. A third parameter, e, is 
an efficiency coefficient for photosynthesis. 

Using this model of leaf growth, we have simulated the 
growth of five tree leaves, which are shown in Fig. [3J to- 
gether with their parameters in the cost function. The 
corresponding thermodynamic functions are shown in 
Fig. |4] The thermodynamic profiles of the tree leaves 
have two regimes that are of particular interest. The 
high energy regime (/3 — >■ 0) can be associated with the 
late stage of the leaf development process, while the low 
energy regime (/3 — ?> cxD)with the early stage of develop- 
ment. 

The internal energy profiles show that, in the early 
stages of growth (/3 — > oo), similar work has been done 
to 'build' all five leaves. In the later stages of develop- 
ment, the amount of work for leaf (d) becomes the lowest, 
which agrees with intuitive expectations, since this leaf 
is the smallest in the sample. The entropy and inter- 
nal energies of leaves (a) and (b) are essentially indistin- 
guishable during the early stages of growth, as are those 
of leaves (c) and (e) . However, discernible differences in 
these quantities appear at later stages. Leaf (d) is dif- 
ferent from all other leaves because of its size, and this 
is also apparent from its thermodynamic functions. The 
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FIG. 4: The internal energy, Helmholtz free energy, and the 
entropy for each of the simulated leaf networks in Fig. |3] 



entropy of leaf (d), in particular, shows a qualitatively 
different profile from all other leaves, which is consistent 
with our earlier interpretation of entropy as the distri- 



bution of energy between the different generations in a 
network. 

The free energy profiles allow us to classify the sample 
leaves into two main groups: {a,b} and {c,d,e}. This 
grouping is confirmed by the low-energy limit (/3 — >■ oo) 
of the entropy profiles of the tree leaves, where the en- 
tropies of the tree leaves correspond to the areas of their 
triangulations in the first generation, as noted in (|5|. 
For leaves with the same initial entropy, differences arise 
when later generations of the triangulation are included. 
This is due to differences in the optimization parameters 
used for the simulations. 

We have described a method for analyzing networks 
based on a triangulation that associates an energy with 
each generation of the network. This enables a partition 
function to be constructed from which standard thermo- 
dynamic functions can be directly computed. The in- 
terpretation of these functions for networks provides a 
powerful classification scheme. The comparison of the 
internal energy, entropy, and Helmholtz free energy for 
several simulated leaves, all of which grow by the same 
mechanism, but with differently parametrized cost func- 
tions, illustrates how our methodology is able to discrim- 
inate between networks that show subtle structural dif- 
ferences. Our method differs from other thermodynamic 
formulations of networks |1HI12J in that our method char- 
acterizes individual networks through their triangulation, 
rather than requiring an ensemble of networks. 

Perhaps the most far-reaching aspect of our analysis 
is that we are able to study the history of a family of 
networks, especially where the early stages of develop- 
ment are similar, but then deviate at some later stage, 
either from an intrinsic mechanism, as shown here, or 
from an external influence. Hence, our method is emi- 
nently suitable for examining biological networks, where 
the entire developmental history, rather than simply the 
fully-developed network, has fundamental importance. 
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lytics LLC. 
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